Analytical solution of the bosonic three-body problem 
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We revisit the problem of three identical bosons in free space, which exhibits a universal hierarchy 
of bound states (Efimov trimers). Modelling a narrow Feshbach resonance within a two-channel de- 
scription, we map the integral equation for the three-body scattering amplitude to a one-dimensional 
Schrodinger-type single-particle equation, where an analytical solution of exponential accuracy is 
obtained. We give exact results for the trimer binding energies, the three-body parameter, the 
threshold to the three-atom continuum, and the recombination rate. 
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The physics of ultracold atoms continues to gener- 
ate immense interest due to the exceptional tunability 
and perfection of the resulting many-body systems 
A particularly interesting topic, first raised in the con- 
text of nuclear physics, concerns three-body bound states 
(trimers) of identical bosons. In seminal work done sev- 
eral decades ago, Efimov 0, H, showed that there ex- 
ists a universal hierarchy of trimer states in the resonant 
case; for a recent review, see Ref. Q. While these states 
were never observed in nuclear physics, first experimen- 
tal reports of Efimov states in the context of cold atoms 
have been recently published (but see Ref. 0] for a 
different interpretation of these results), leading to re- 
newed interest in Efimov physics also among theorists 

Q. For instance, when atoms are confined in a trap, 
confinement effects can affect Efimov states in interest- 
ing ways [13, [m, ■ More generally, Efimov physics has 
strong experimental relevance since it imposes limits on 
the atomic cloud lifetime via three-body recombination 
processes. 

Here we reconsider the homogeneous three-dimensional 
(3D) case for three identical bosons and present a novel 
solution of the Efimov problem, which is simpler than 
existing approaches 0] and thus allows for analytical 
progress. Our scheme is based on a mapping of the three- 
body integral equation to an effective ID Schrodinger- 
type single particle equation, where quantum-mechanical 
intuition is available. Zero-range interactions for three 
bosons is an ill-defined problem, with a dense [isl. [T3 | and 
unbounded spectrum ('Thomas collapse') [l^ of bound 
states. In addition to the usual two-body scattering 
length a [l[ , the regularization of the three-boson prob- 
lem thus requires an additional parameter R* > related 
to the finite range of the interaction potential. We here 
consider the case of a narrow Feshbach resonance, where 
the short-distance regularization does not involve higher 
angular momentum {I > 0) partial waves 0]. Remark- 
ably, the resulting Schrodinger-type equation for the ID 
motion of a fictitious quantum particle in a certain po- 
tential, see Eq. pd|) below, can be solved analytically. 
As is shown below, corrections to this solution are ex- 
ponentially small at low binding energies. Our theory 



reproduces the known binding energy hierarchy for shal- 
low trimer states, sec Eq. ([Ti)) . We provide exact re- 
sults for three key quantities: (i) the three-body parame- 
ter K*, see Eq. p^ . which determines many observables 
known to be universal functions of Q; (ii) the scat- 
tering length a'^ where the Efimov states disappear into 
the three-atom continuum, see Eq. ([19]); and (iii) the re- 
combination rate ccrec for three-body recombination to 
a weakly bound level, see Eq. (|2ip . Our result for 
depends on short-distance properties and is specific to 
narrow resonances, while the a'^ and arec results apply 
for any type of resonance. Finally, our approach should 
also simplify future calculations in confined geometries. 

To begin with, let us derive a suitable three-body 
integral equation encoding Efimov physics. Inspired 
by Ref. [16|, we use a two-channel model to regular- 
ize the theory, where the boson operator Ok, describ- 
ing a momentum-k atomic state, is coupled to the bo- 
son operator 6k describing a closed channel (molecular) 
state of momentum K and energy offset Eq. This model 
reproduces basic features of Feshbach resonances 
llSf . which are commonly used to tune the interaction 
strength in experiments. With coupling strength A 
and efc — h'^k'^/2m, the second-quantized many-particle 
Hamiltonian for a narrow Feshbach resonance is 

H = ^e;caUk + 5^(i?o + ei^/2)6j,6K (1) 
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The two-body problem is solved by the Ansatz = 
(/^^K=o + Sk ^kakaLk)|0), where the Schrodingcr equa- 
tion {H ~ E)\'ip) = imphes coupled equations for the 
complex coefficients /? and A^. [isj . 



k 

AP + {2ek - E)Ai, = 0. 



(2) 



The second equation is solved by an incoming scattering 
state of momentum ko and energy E = 2eko ~ h^k^/m, 
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where = {2n)^d{k - ko) + 

body scattering amplitude / = 
fine 

1 1 2TTh'^EQ 



with the two- 
-mA/3/(47rft2). We de- 



R* 



(3) 



a £ toA^ 
where £ is a short-distance cutoff len gth [l9| of the order 
of the van der Waals potential size |16|. For a narrow 
Feslibach resonance, one has R* ^ i [8|,[20[. Equation ^ 
results in /(fco) = —[a~-^+iko+R*kQ]~^, which coincides 
with the familiar effective range expansion [l[ , but is not 
restricted to the regime |fcoi?*| <C 1. The three-body 
problem can then be solved using the Ansatz 
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without the need for any additional rcgularization. The 
(unsymmetrized) three-boson wavefunction Ak.Ic can be 
taken as even function of k. The Schrodinger equation 
again leads to coupled algebraic equations, 

Eo-E + ^^)Pj^ (4) 
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2^k-K/2,-k/2-3K/4) — 0: 
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Akm = 0. 



Since we search for bound states, we now put E = 
—h^X^/m < 0. The second equation is then solved by 
^K,k = -imAPK/h'^)[3K'^ /4: + P + \^]-\ while the first 
determines the three-body scattering amplitude /3k- Em- 
ploying the definitions of a and R* in Eq. we find 



{3R*Ky4 + Lx - l/a)pK = 0, 
where we introduce the operator 



(5) 
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/\ 2 + K' • K + A2 ■ 



This is the three-boson integral equation, cp. Refs. 
Hi [IH- In contrast with the commonly used Bethe- 
Peierls boundary conditions or pseudo-potentials [2l| , the 
Hamiltonian ([1]) is manifestly self-adjoint, and thus our 
two-channel description provides a natural regularization 
scheme free of the mathematical difficulties encountered 
otherwise [H, • It differs from the real-space scheme 
used originally by Efimov 0, H, Q , but leads to the same 
universal trimer energy hierarchy. Rotational symmetry 
expresses the solution of Eq. ([5]) in terms of a 'wave- 
function' il]{k) for a fictitious ID problem, /3k = i'{k)/k, 
with k = |K| > 0. We thus integrate over the angular 
variables in Eq. ([5]), and obtain 

(v/3fcV4TA2 -a-i +i?*(3fcV4 + A2)) 7/.(fc) (7) 
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^(fc') = 0. 



Since the logarithmic term is odd under k' —k', we 
extend k to include negative values, but require = 
—il;{k). The integral in Eq. ^ is then taken over all k', 
with the replacement 2/t: I/tt. 

We next show how Eq. ([7]) can be mapped to a 
quantum-mechanical single-particle equation. To that 
end, we employ the transformation 



2A 
7! 



sinh^, ^ e 



(8) 



and search for antisymmetric solutions "0(0 — ^i-'i^O- 
Substitution ([5]) is helpful because (i) the square root 
in Eq. ([7]) rationalizes, and (ii) the logarithmic kernel 
becomes homogeneous. To sec that, we decompose the 

logarithmic factor in Eq. ([7]) as In (^fsf^fj^^jT^pj^^ + 

In ^ ^2(;+;')^ee+e'+i ) ■ I^'^^b terms yield the same contri- 
bution in Eq. ([7]), which is seen by letting ^' — ^' in the 
integral associated with the second term and exploiting 
antisymmetry, = —ipiC)- This leads us to the 

difference kernel T(^) 
with Fourier transform 
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An 8 sinh(7rs/6) 
3V3 s cosh(7rs/2) ' 



(9) 



This kernel acts on a test function g{^) as a differen- 
tial operator, J^^d^Ti^ - ^gi^) = T (-i^) giO- 
The function ([9]) thus plays the role of a kinetic energy 
operator. For the standard non-relativistic Schrodinger 
equation, T{p) = p^/2m; our notation s for 'momentum' 
follows convention [Sj. Here the dispersion relation ^ 
starts as T(s) ^ at small momentum and levels off to 
47r/ (3v^) as s ^ c». It is thus bounded from below and 
from above, similar to what happens for the typical band 
structure of a solid. Note that we slightly abuse terminol- 
ogy here, since ^ parametrizes the physical momentum, 
see Eq. ([S]), and hence s in reality corresponds to a spa- 
tial variable. Nevertheless, this analogy is quite helpful 
to make contact with elementary textbook quantum me- 
chanics, and we will henceforth denote s as 'momentum' 
conjugate to the ID 'space' variable ^. We also introduce 
the symmetric single-particle potential 



1 1 

aX cosh^ 



i?*A cosh^. 



(10) 



and the 'energy' £ ^ 47r/(3V3) - 1 ~ 1.41839. After 
rescaling cosh(^), Eq. ([7]) then assumes 

the equivalent final form 



T[^^-]+UiO-£ 



m = 0, 



(11) 



subject to the antisymmetry condition '(/'('?) = ~0(^O- 
This equation formally describes the ID motion of a 
quantum particle with non-standard dispersion relation 
^ in the potential U{£,) at energy £. 
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Let us now consider the resonant limit 1/a = 0. The 
solution of Eq. ([TT]) for iZ* = is given by V(0 ^ 
sin(soC)i where sq ~ 1.00624 is the real and positive root 
of T(so) = £■ However, this does not lead to well-defined 
energy eigenstates d, [H, [l^ . We thus need to analyze 
the regularization mechanism due to i?* > 0, which is 
quite simple to understand within our picture: Potential 
pn)) approaches +oo at ^ ^ ±oo, and hence all eigen- 
states must be quantized bound state solutions, similar 
to what happens for a simple harmonic oscillator. In our 
case, the 'energy' £ is always fixed but the true spectral 
parameter is A: only those values of A are allowed (possi- 
bly a finite set or countable infinity), where a bound state 
with energy £ exists. These discrete values A„ (integer 
n) then determine the Efimov trimer bound state ener- 
gies En = —h'^X'^/m. For R*X <C 1, we have n » 1, and 
zero energy represents a spectral accumulation point, see 
Eq. dm) below. Taking ^ > with = --0(0, the 

potential U{£,) can be neglected against £ in the region 
C < where = ln[2/(i?*A)] > 1. Here the solu- 
tion must therefore be = ci sin(soO, with complex 
amplitude ci . On the other hand, for all ^ ^ 1 (which in- 
cludes ^ ~ ^*); the potential takes the form U{£,) — e^~?* . 
Shifting ^ by this regime is thus described by a uni- 
versal (parameter-free) equation 



T 



m = 0. 



(12) 



Note that antisymmetry of V'(C) should not be imposed 
here because of the shift. For ^ ^ — oo, the exponen- 
tial term in this equation vanishes, and the asymptotic 
behavior ip{(,) sin(soC + ""7) with a non-trivial phase 
shift 7 is expected. Coming back to the original ^, the 
solution for 1 < ^ < ^, is -02 (■?) = C2 sin[so (C - ) + "'7] , 
where C2 is another amplitude. Evidently, both 'ip2 and 
ipi should match within the broad region 1 ^ ^ <C ■f*. 
This implies the quantization condition 



e,(A„) = ln[2/(i?*A„)] = Tr{n + j)/so, 



(13) 



yielding the on-resonancc ERmov trimer energies 

E„ = te-2™/.o (14) 

TO 

with the famous universal ratio En+i/En = e^'^'^l^^ ~ 
1/515.03 between subsequent levels. Note that i?*A ~ 



is exponentially small for shallow Efimov trimcrs. 



Since all corrections to our approach arc of order i?*A, 
deviations from our results (e.g., |(5£'„)/i?„|) are expo- 
nentially small for n — s- oo. The three-body parameter 
K* 0, defined only up to multiplicative factors of e^/^°, 
is here taken as k^R* = 2e~'"'^^'^°. To obtain it, we now 
calculate 7 from the universal problem (|12[) . 

Remarkably, the exact solution of Eq. (|T^ can be ob- 
tained as a Barnes-type contour integral [22 HI] , 



ioo+0+ 



ioo+0+ 



(15) 



where Eq. (fT2|) implies the recurrence relation [T{iv) — 
£]C{iy) = -C{v + 1). Applying f - £ to Eq. ^ 
thus changes the integrand to —e~'^^C{v + 1). The de- 
sired result — 6^-0(0 is obtained if no pole is crossed 
when the contour is shifted back, v ^ v — \, see also 
Ref. 0. Therefore C{v) must have no poles within 
the strip < Yle{v) < 1 of the complex plane. Us- 
ing the Weierstrass theorem [22|, we can now express 

T{iu) — £ = np=o v'^-b^ ^ ^ convergent infinite prod- 
uct in terms of its poles ±6p, where bp = 2p + 1, and 
its zeroes ±Up. Two zeroes are on the imaginary axis, 
uo ~ iso, while all others are real and correspond to 
ui =4,-^2 = 4.6, ... [1, 0]. The function satisfying the 
recurrence relation and the required analytic properties 
is C{l') = [tt/ sin(7r(i/ — iso))] C+{i'), where 



c+(^)=n 



T{iy 



bp) 



p=o + -^ + Up) 



(16) 



with the gamma function T{z). In fact, C(i/) has no 
poles in the strip < Rc(i') < 2, while the poles at 
i> ^ 2 and v ~ ±iso imply from Eq. ()15|1 the asymptotic 
behaviors V'(C) ^ e~^^ for ^ — > 00 and, as expected, 
"0(0 ^ sin(so'^ + T^"f) for ^ — > —00. The phase 7 now 
follows from the ratio of the residues at the two poles 
1/ = ztisQ. We find the exact result 



7 = ^ - -ArgC+(iso) 

Z TT 



-0.090518155. 



Hence the three-body parameter is given by 
K^R* = 2e"''''/"« ~ 2.6531. 



(17) 



(18) 



This exact result roughly agrees with the available nu- 
merical estimate w 2.5 0, We have also per- 
formed a more accurate numerical check [23|- This gives 
K^R* ~ 2.653(1), in good agreement with Eq. (fT8|) . 

Using Eq. (fT7|) we can also compute the scatter- 
ing length a'^ < where the Efimov state joins the 
three-atom threshold at A ^ 0. In that limit, the 
off-resonant potential p0|) can be written as U {£,) = 
(2/|a|A)e~^ + (i?*A/2)e^, and the same phase shift 7 thus 
determines the behavior at ^ ^ ±00. Equation p3p 
then implies ln(2/i?*A) ^ 7r(?i + -/)/sq and ln(2/|a|A) = 
7r(7i' — 7)/so (with integers n,n'). This in turn gives 
a'^/R* = -e^'^'^/'^f (modulo factors of e''/^°), and hence 
the exact result for the three-atom continuum threshold. 



T7/S0 



-1.50763, 



(19) 



consistent with the previously obtained numerical result 

= -1.56(5) g. 

Finally, let us briefly address the three-body recom- 
bination rate aj-ec- For energy E = Q'^ and a > 0, 
the solution of Eq. @ proceeds as before, but now in- 
cludes a source term in Eq. ([5]) . Fourier transforma- 
tion to real space and rescaling r \/3r/2 as well as 
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/3(r) 6V27ri?*a(l + 4/3(r)/\/3) results in 

(-i?*Ar4--|z,„i)^(r).i, (20) 

where Lq is given in Eq. ([6]). For r ^ oo, the so- 
lution of Eq. ([20]) contains an outgoing spherical wave 
Ae'^'^^ /r^ describing atom and dimcr separating af- 
ter recombination, where h^k'^/m is the weakly bound 
dimer energy calculated from f~^(ikd) = 0. Follow- 
ing Ref. [l6[, Orec is then obtained in the form a^-^c = 
96V37r2(fioVTO)|A|2(l -f 2R*kd). For R* > a, Lq be- 
comes negligible in Eq. ([20]) . and the result aroc = 
192V3Tr'^{h/m)Va7B* follows. For R* < a, on the 
other hand, a complete analytical solution is also possi- 
ble by relating Eq. ([20|) to the universal problem (fT2|) . 
In fact, in the broad region R* r a, Eq. ([20|) is 
basically the Fourier transformed equation p2[) . leading 
to /3(r) ~ sin[so ln(r/i?o)]j with the renormalization 
length 3 i?o = R*eM'r+im+Arer{iso)]/s„ ^ o.55890i?*. 
This imposes a boundary condition for R* <Si r <Si a on 
solutions to Eq. (|20)) . which in turn allows to compute 



A and hence arec in closed form. Postponing details to 
elsewhere, the final result is, cf. Refs. [g. [25j. 

2^71^(4^-3^3) ha'' . 2 

arec = —^7 ^ sm dr (21) 

smh (ttso) 

with = (5i - Arg(l + e-2''"«e2^'^i): where (5i = 
So ln(aK*/2) + tt/2 — nj. 

To conclude, we have presented a fresh theoretical ap- 
proach to the Efimov problem of three identical bosons. 
It is based on mapping the three-body integral equation 
for the scattering amplitude to a Schrodinger-like equa- 
tion for a quantum-mechanical particle in a ID potential. 
For low binding energies, this permits a solution of ex- 
ponential accuracy, and we have derived exact results for 
several key observables. 
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